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We present the first theoretical treatment of the formation of highly excited neutral H atoms (H*) 
for strongly driven H2. This process, with one electron ionizing and one captured in a Rydberg state, 
was recently reported in an experimental study in Phys. Rev. Lett 102, 113002 (2009). We show 
that two pathways underlie this process: one mostly of non-sequential nature — more so for small 
intensities — resembling non-sequential double ionization, and one mostly of sequential nature — more 
so for high intensities — resembling double ionization through enhanced ionization. We also predict a 
new feature, asymmetric energy sharing between H* and H + with increasing intensity. This feature 
is a striking demonstration of the influence the electron has on the nuclear motion. 

PACS numbers: 33.80.Rv, 34.80.Gs, 42.50.Hz 



A wealth of physical phenomena is manifested dur- 
ing fragmentation of molecules driven by intense infrared 
laser fields. In the case of H2 — the most fundamental di- 
atomic molecule — such phenomena include, for instance, 
non- sequential double ionization (NSDI) [TJ [2] with the 
re-collision (three-step) model [3] underlying this process; 
double ionization through enhanced ionization (EIDI) 
with one electron escaping early on while the other ion- 
izes later at intermediate distances of the nuclei of the 
remaining molecular ion [4-6 . Very recently, formation 
of highly excited neutral H atoms [7 was identified as 
yet another interesting phenomenon in driven H2. It was 
attributed to the "frustrated tunnel ionization" mecha- 
nism leading to the recapture of an electron. First seen 
in atoms [8 this mechanism appears to be a general 
one underlying the formation of excited fragments dur- 
ing Coulomb explosion in strongly driven molecules [9]. 
However, a full theoretical treatment and understand- 
ing is pending. This is partially due to the fact that 
theoretical studies including both nuclear and electronic 
motion in strongly driven systems are computationally 
very demanding. Most studies fix the nuclei and focus 
on electronic motion [TJ [10] or ignore the electronic con- 
tinuum and study nuclear motion [11] ; in a recent refined 
study the electronic structure is used as input to study 
the energy spectra of the Coulomb exploding nuclei [12] . 

In this Letter, we present the first comprehensive the- 
oretical study on the break-up of strongly driven H2. We 
use a three-dimensional semi-classical approach that ex- 
plicitly accounts for the Coulomb singularity. We focus 
on the Coulomb explosion process involving a highly ex- 
cited neutral hydrogen atom (H*) and a proton. The dif- 
ficulty in its theoretical treatment lies in the fact that, as 
we show, two-electron effects pervade the two pathways 
underlying the formation of H*. This is in contrast to 



one-electron effects pervading the well known enhanced 
ionization process. Moreover, correctly accounting for 
this process requires treating the electronic and nuclear 
motion on an equal footing. By doing so, we identify 
a transition with increasing laser intensity from sym- 
metric to asymmetric energy sharing between H + and 
H*. We show this asymmetry in energy sharing to be a 
striking signature of how the motion of a light particle — 
electron — influences the motion of heavy ones — nuclei. 
This questions whether this asymmetry could be pre- 
dicted by the widely used in quantum mechanics Born- 
Oppenheimer (BO) approximation, since the latter as- 
sumes that the massive nuclei are nearly fixed with re- 
spect to electron motion. 

Model. Our calculations are performed in a classical 
framework with the Coulomb singularity explicitly ac- 
counted for. Previously, we successfully used a three- 
dimensional classical model for fixed nuclei diatomics to 
predict new phenomena: an antiparallel electron escape 
for intermediate laser field intensities in the over-the- 
barrier regime [13]. The current model extends our pre- 
vious technique to also account for nuclear motion. 

Our model entails setting up the initial phase space dis- 
tribution. The electronic initial state is detailed in [13] 
both for the tunneling and the over-the-barrier regime. 
We use 0.57 a.u. and 1.28 a.u. as the first and sec- 
ond ionization potentials, respectively. We take the ini- 
tial vibrational state of the nuclei to be the ground 
state, Eq w 0.01 a.u, of the Morse potential Vm(R) — 
D(l - e -^ R - R ^) 2 with D = 0.174 a.u., f3 = 1.029 
a.u. 14] and R the inter-nuclear distance. We choose 
the Wigner distribution to describe the initial state of 
the nuclei due to the good agreement with the position 
and momentum distributions of the initial quantum me- 
chanical state. (For the Wigner function of the ground 
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state Morse potential see [II]). We consider an intensity 
in the tunneling regime of 1.5xl0 14 Watts/cm 2 and an 
intermediate intensity in the over-the-barrier regime of 
2.5xl0 14 Watts/cm 2 . For this latter intensity, for the 
majority of trajectories the electronic state is initiated 
with the tunneling model, thus the term intermediate. 
Both intensities are high enough to justify restricting the 
initial distance of the nuclei to the equilibrium distance 
of the ground vibrational state, Rq = 1.4 a.u. [15]. In 
our study, the laser pulse is aligned with the molecular 
axis, has frequency uj = 0.057 a.u. and a 10 cycle dura- 
tion plus a 2 cycle turn-off — we use a cos pulse with a 
cos 2 turn-off. After setting up the initial state, we trans- 
form to a new system of coordinates, the so called "regu- 
larized" coordinates [16] . This transformation explicitly 
eliminates the Coulomb singularity [T3l [17] . We regular- 
ize using the global regularization scheme described in 
[15] . Finally, we propagate in time using the Classical 
Trajectory Monte Carlo method [19]. (We note that to 
correctly account for EIDI and formation of H* we al- 
low for both electrons to tunnel during propagation in 
a manner similar to that outlined in [20].) We now se- 
lect those trajectories leading to a break up of H 2 with 
H + , H* (where * denotes an electron in a n > 1 quan- 
tum state) and a free electron as fragments. To identify 
the electrons captured in a Rydberg n state we first find 
the classical n r 



l/y / 2\E\ where E is the potential plus 
kinetic energy of the electron. Next, we assign to this 
classical number a quantum one such that the relation- 
ship ((n-l)(n-l/2)n)^ 3 <n c < (n(n^l/2)(n^l))^ 3 , 
derived in [21], is satisfied. 

As a result of our calculations we find that this H* 
ionization channel in the strongly driven H2 accounts for 
roughly 7-10% of the break-up events. We compute the 
probability distribution as a function of the kinetic en- 
ergy of either one of the fragments H + , H* for 1.5 xlO 14 
Watts/cm 2 and 2.5xl0 14 Watts/cm 2 , see Fig. [l] a) and 
b). We find that the maximum of the probability distri- 
bution is around 3.5 eV. This is in very good agreement 
with the experimental results in [7 and suggests that 
our model accurately describes this process. (The ex- 
perimental peak at about 0.5 eV due to bond-softening 
is not addressed by the current work). Furthermore, we 
find that the distribution of principal quantum numbers 
n in H*, peaks around n = 8 resembling results for the 
atomic case [8]. 

In a first attempt towards understanding, in [7\ the for- 
mation of H* was attributed to one electron effects, i.e., 
"frustrated" (tunnel) ionization of the electron in — 
assuming the first electron ionizes immediately. "Frus- 
trated" ionization refers to the process where after tun- 
neling an electron quivers in the laser field; when the 
field is turned off the electron does not have enough drift 
energy to overcome the Coulomb potential and gets cap- 
tured in a Rydberg state. Our full classical calculation 
shows this conjecture to be partly true. It reveals that 
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FIG. 1. (Color online) Top row: Probability distribution for 
either of H + or H* to have a certain amount of energy in the 
H* formation channel; Bottom row: The double differential 
probability in energy for H + and H*. Left column at 1.5x 10 14 
Watts/cm 2 and right column at 2.5xl0 14 Watts/cm 2 . 



two electron effects pervade the two pathways that we 
find underlying the formation of H* . In the pathway we 
refer to as first, while the first electron — the one that 
tunnels in the initial state — quivers in the laser field, the 
initially bound second electron gains energy and ionizes. 
The second electron can gain energy through two mecha- 
nisms a) a non-sequential one, where the first tunnel elec- 
tron transfers energy through a re-collision; b) a sequen- 
tial one, where the the initially bound electron, indepen- 
dently of the first tunnel electron, ionizes at an interme- 
diate distance of the nuclei through enhanced ionization 
of H^. Both for a) and b) the first electron eventually un- 
dergoes "frustrated" ionization. In both cases the timing 
of the ionization of the second electron determines the 
kinetic energy of the fragments. 

In the second pathway, the initially bound second elec- 
tron gains energy through two mechanisms a) a non- 
sequential one where it undergoes a re-collision with the 
quasi-free first electron; b) a sequential one where the 
first electron ionizes fast and the bound electron gains 
energy through the lowered potential due to Coulomb 
explosion of the nuclei. Both for a) and b) the initially 
bound electron undergoes "frustrated" ionization. 

Two-electron effects mostly prevail in the first path- 
way. Indeed, in this pathway the first electron quivers 
in the proximity of the nuclei before the bound elec- 
tron ionizes. However, in the second pathway the first 
electron ionizes early on while excitation of the bound 
electron follows only afterwards. A finding corroborat- 
ing the significant role of two-electron effects in the first 
pathway for small intensities is that the contribution 
of the first pathway decreases from 7% of all break-up 
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events for 1.5 xlO 14 Watts/cm 2 to half that for 2.5 xlO 14 
Watts/cm 2 . This is consistent with the diminishing role 
of re-collisions with increasing intensity. On the other 
hand the second pathway's contribution remains con- 
stant around 4% for both intensities. Moreover, using 
as a rough criterion the presence of a re-collision (mini- 
mum in the distance of the two electrons) we can further 
separate the trajectories of the first and second path- 
ways to sequential and non-sequential ones. We find that 
the non- sequential/sequential trajectories correspond to 
the first electron tunneling at positive/negative phases of 
the laser field. While non-sequential trajectories prevail 
for 1.5 xlO 14 Watts/cm 2 sequential ones prevail for 2.5 
xlO 14 Watts/cm 2 . Typical trajectories of the first and 
second pathway are shown in Fig. [2] b) and Fig. |3]b), 
respectively. 

From the above we conjecture that the non-sequential 
part of the first pathway of H* formation for 1.5 xlO 14 
Watts/cm 2 is the one that resembles the most non- 
sequential double ionization (NSDI) — for H* double ion- 
ization is frustrated and the first electron remains cap- 
tured. Indeed, comparing Fig. [2] a) with Fig. [| c) we 
find that the probability distribution for the energy of 
H + or H* in the non- sequential part of the first pathway 
of H* is similar to the probability distribution for the 
energy of the two H + fragments in the Delayed pathway 
of NSDI. We find that NSDI is roughly 2% of all break- 
up events. As is well known, the Delayed pathway — one 
electron ionizes with a delay, more than a quarter laser 
cycle, after re-collision — is one of the major energy trans- 
fer pathways in NSDI. (For details on the separation of 
NSDI pathways see [13].) A comparison between Fig. [2] 
a) and Fig.[2]c) further shows that in the first pathway of 
H* formation the bound electron ionizes later compared 
to its ionization time in the Delayed pathway. Indeed, 
larger ionization times correspond to the nuclei Coulomb 
exploding at a later time, resulting in final energies of 
the nuclei smaller for H* compared to NSDI. 

The sequential part of the second pathway of H* for- 
mation for 2.5 xlO 14 Watts/cm 2 is the one resembling 
the most double ionization through enhanced ionization 
(EIDI) — for H* double ionization is frustrated and the 
second electron remains captured. As is well known (see 
[2]) in the EIDI pathway the first electron ionizes fast 
while the initially bound one ionizes at an intermediate 
distance of the nuclei. We find that the EIDI process 
accounts for more than 50% of all events in agreement 
with previous results [2]. For EIDI the probability dis- 
tribution for the energy of H + shown in Fig. [3] c) peaks 
at 3.5 eV, also roughly in agreement with previous re- 
sults, see for example [12]. A comparison of Fig. [3] a) 
with Fig. [3] c) shows that indeed the sequential part of 
the second pathway of H* formation resembles EIDI. 

We have now established the pathways of energy trans- 
fer underlying H* formation. Next, we seek further in- 
sight into the underlying processes by considering the 
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FIG. 2. (Color online) For 1.5 xlO 14 Watts/cm 2 a) Probabil- 
ity distribution for either H + or H* to have a certain amount 
of energy in the non-sequential part of the first pathway of 
H* ; c) Probability distribution for either of the H + fragments 
to have a certain amount of energy in the NSDI process: cir- 
cles show the total probability distribution while the triangles 
show the contribution from the Delayed pathway; b) A typical 
trajectory for the non-sequential first pathway of H* forma- 
tion with the solid/dashed black line depicting the nucleus 
escaping in the same/opposite direction of initial tunneling, 
the dashed green line the bound electron and the red line the 
first electron. In the inset of Fig. [2]b) we show how the first 
electron undergoes "frustrated" ionization. 
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FIG. 3. (Color online) For 2.5 xlO 14 Watts/cm 2 a) same as 
Fig. p a) but for the sequential part of the second pathway 
of H formation; c) Probability distribution for either of the 
H + fragments to have a certain amount of energy in the EIDI 
process; b) same as Fig. [5] b) but for the sequential part of 
the second pathway. In the inset of Fig. [3] b) we show how 
the second electron undergoes "frustrated" ionization. 



energy sharing between H + and H*, see double energy 
differentials in Fig. [l]c) and d). We find that H + and H* 
share the energy equally for 1.5xl0 14 Watts/cm 2 and 
unequally for 2.5xl0 14 Watts/cm 2 , albeit the difference 
in energy sharing is small. The equal energy sharing 
is in very good agreement with the experimental results 
in [Tj for an intensity of 3xl0 14 Watts/cm 2 . However, 
the asymmetry we predict in energy sharing at 2.5 xlO 14 
Watts/cm 2 suggests that the experimental results have 
contributions from smaller intensities. 

We find the transition from symmetric to asymmetric 
energy sharing to be related to the phase of the field 
at the time the first electron tunnels. For 1.5 xlO 14 
Watts/cm 2 the distribution of the phase of the laser field 
when the first electron tunnels is centered around zero, 
see Fig. [4] c). According to the three-step model, for 
phases around 0° the first electron returns close to the 
nuclei. Indeed, as shown in Fig. [4] a), the first electron 
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oscillates between the two nuclei in the z (field) direction 
for several periods before it finally gets captured in a 
Rydberg state of H. The return of the first electron close 
to the nuclei favors transfer of energy through the non- 
sequential H* mechanism, indeed the case for 1.5 xlO 14 
Watts/cm 2 as discussed above. Moreover, with the first 
tunnel electron oscillating between the two nuclei, the 
nuclei experience the same attractive force from the elec- 
tron. This force combined with the repulsive force due to 
Coulomb explosion, results in a symmetric energy shar- 
ing of H + and H*, see Fig. [l]c). 




field phase (degrees) 

FIG. 4. (Color inline) Top row: For the first pathway, mean 
value of the position component z of nuclei and electrons when 
z > for the second electron; lines are the same as for Fig. [2] 
b). Bottom row: Probability distribution for both pathways 
for the phase of the laser field when the first electron tunnels. 
Left/right column for 1.5xl0 14 /2.5xl0 14 Watts/cm 2 . 

For 2.5 xlO 14 Watts/cm 2 the distribution of the phase 
of the laser field when the first electron tunnels is cen- 
tered mostly around -30°, see Fig. gd). For large neg- 
ative phases the first electron does not return close to 
the nuclei. Instead, it remains on the same side that it 
tunneled from in the z direction before it gets captured, 
see Fig. |4]b). This behavior does not favor re-collisions 
and is consistent with our finding of sequential trajecto- 
ries prevailing the formation of H* at this higher inten- 
sity. Moreover, the behavior of the first electron causes 
an asymmetry in the force experienced by the Coulomb 
exploding nuclei. The nucleus escaping along the same 
direction as the first electron experiences from this elec- 
tron an attractive force that adds to the force due to 
Coulomb explosion. However, the nucleus escaping in 
the opposite z direction experiences an attractive force 
from the first tunnel electron that decreases the force 
from Coulomb explosion. This explains the asymmetry 
in energy sharing for 2.5 xlO 14 Watts/cm 2 in Fig. [l]d). 

In conclusion, we have shown that the H* formation 
in the break-up or Coulomb explosion of H2 in strong 
laser fields proceeds through two pathways. One where 



the first electron undergoes "frustrated" ionization while 
the initially bound electron ionizes by absorbing energy 
from the field or via a collision process. In the other 
pathway, the first electron ionizes and the second electron 
gains energy either from the field or the first electron 
and undergoes "frustrated" ionization. Furthermore, we 
have shown that for higher intensities electron tunneling 
at large negative phases gives rise to asymmetric energy 
sharing between H + and H*. This asymmetry might not 
be possible to predict using the BO approximation since 
the latter assumes that the nuclei are nearly fixed with 
respect to the electronic motion. This could be checked 
quantum mechanically with since the asymmetry in 
energy sharing is mostly a one electron effect. 
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